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Abstract 

In the present work, we implement an explicit two-loop renormalization 
of a two-dimensional flat Fermi surface (FS) in the framework of a field 
theoretical renormalization group (RG) approach. In our scheme, we derive 
the RG equations for both coupling functions and Fermi energy. In this way, 
we are able to probe the existence of spin-charge separation by showing that 
the low-energy sector of the system is in fact a non-Fermi liquid. In addition, 
associating the true interacting FS to the infrared stable (IR) fixed point 
of the Fermi energy, we demonstrate here that it either acquires a small 
curvature and behaves as a "Luttinger liquid" or it suffers a truncation in 
k-space depicting an insulating spin liquid. 

1 - Introduction 

A better understanding of the physical properties of highly interacting electrons 
in two spatial dimensions (2d) is central for high-Tc superconductivity. Soon 
after the discovery of the high-Tc superconductors, Anderson [I] suggested that a 
strongly interacting 2d electron gas should resemble a Id Luttinger liquid state. 
This question remains unresolved to this date. Thanks to the high precision of the 
angular resolved photoemission experiments performed in a variety of materials [2 , 
we know, at present, important facts concerning the Fermi surface (FS) of the 
cuprates. The FS's for underdoped and optimally doped Bi2212 and YBaCuO 
compounds contain both flat and curved sectors j^j. As a result, they are nearly 
perfectly nested along certain k-directions. As is well-known, whenever there is a 
flat FS, the corresponding one-electron dispersion is ld-like in momentum space. 

Originally, the cuprates are Mott insulators which become metallic at very low 
doping At half-filling, Hubbard-like models have a square shape FS imposed 
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Figure 1: The initial Fermi surface. The corners are rounded to avoid van Hove singu- 
larities. 



by electron-hole symmetry. The FS changes as we vary the filling factor and, as 
soon as it is lightly doped, it acquires nonzero curvature sectors in k-space. In the 
immediate vicinity of half-filling, there are, at most, isolated curved spots in mo- 
mentum space. Consequently, in a zeroth order approximation, one may neglect 
their presence altogether. Following this scheme, several workers investigated the 
properties of a 2d electron gas in the presence of a totally flat FS |SJ El E] • In 
their approaches, the FS is always kept fixed and never deviates from its original 
flat form. Besides that, their results conflict with each other. Conventional per- 
turbation theory calculations [3], parquet method results |6., as well as one- loop 
perturbative RG calculations [7j indicate that, for repulsive interactions, there is 
never a Luttinger liquid state in 2d. In contrast, applying bosonization methods, 
Luther was able to map the square FS onto two sets of perpendicular chains jHj • 
As a result of that, the corresponding electron correlation functions become sums 
of power law terms with exponents only differing in form from those of a Luttinger 
liquid [9 . We revisit this problem in this letter. 

We report a two-loop field theoretical RG calculation for the electron gas in 
the presence of the same FS model as used by Dzyaloshinskii and co-workers jHJ- 
The novel aspect of our work, apart from taking into account important higher 
order corrections, is the fact that we show explicitly how the FS changes its shape 
due to interactions. As a result, we are able to determine when the FS may suffer 
a truncation in k-space. To our knowledge, this was never done before in such a 
systematic and detailed form. Needless to say, to show how the FS is renormalized 
by interaction is a very intricate many-body problem [TU] . 

To deal with this here, we calculate how the charge renormalization functions 
Z and all physical parameter vary along the renormalized FS itself, by means of 
appropriate RG flow equations for the Fermi energy and coupling functions. We 
observe that Z is nullified at FS and, as a result, there exists spin-charge separation 
in 2d. In addition, we explore the existence of nontrivial fixed points which vary 
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continuously along the Fermi surface. Accordingly, we show that the renormalized 
FS naturally develops a nonzero curvature, if one associates its physical nature 
to the infrared (IR) stable fixed points. It follows from this that its behavior is 
regulated by the variation of the anomalous dimension exponent 7* with respect to 
the momentum pn along FS. When 1/2 < 7* < 1 the renormalized FS is truncated 
and develops a charge pseudogap in this region of k-space. The physical system 
behaves as an insulating spin liquid in this case. In contrast, when < 7* < 1/2 
the non-Fermi liquid is metallic and resembles a Luttinger liquid state. 

2 - Renormalized hamiltonian and electron self- 
energy 

Our starting point is a strongly interacting 2d electron gas in the presence of the flat 
FS shown schematically in fig. d In order to keep a closer contact with well-known 
works in one-dimensional physics ^T] , we split the FS in four patches. However, we 
expand the bare single-particle energy dispersion in the vicinity of the renormalized 
(i.e physical) FS. The parametrization of the corresponding interactions is shown 
schematically in fig. 2. If we use this model to calculate physical quantities using 
a naive perturbation theory, we find divergent results for particular values of the 
external momenta ^21 EH] • We circumvent this problem following the standard 
field theory procedure of introducing appropriate counterterms in the hamiltonian 
to render the physical parameters finite in all scattering channels jHj . In this way, 
the original hamiltonian is rewritten in a more convenient form in terms of the 
low-energy parameters, which are in turn physically measurable. Thus, we have 
H = Hr + He, where 

Hr = ^v FR (\p ± \ - k FR )ip Ra (p)^R a (p) 

per 

+ + U ™ + U ™ + U ^L(P + q - k)Vi- CT (k)^,-a(q)VMP) (!) 

pqker 

and 

He = £[1^(4* - k lB) + (^(P)— - 1> FR (\P±\ ~ k FR )A(p)^(p) 

— ZTflR TflR 

per 

pqker a 

(2) 
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Figure 2: The parametrization of the interactions in the present model: (a) U±r- 
processes, (b) [^-processes, (c) [/^-processes, and (d) [/^-processes. 

Here the subscripts "R" and "B" stand for renormalized and bare respectively. 
Besides, Z(p) is the charge renormalization function, which is well-defined only 
in the following regions of k-space depicted in fig. 1: —A < p± =F k F R < A and 
—A < p» < A. Finally, the renormalized and bare fields are related to each 
other by iPb(p) — Z 1 ^ 2 (p)4'r(p), whereas the renormalized and bare couplings are 
connected by Yii=i Z^ 2 (Pi)UiB = UiR + AUffi . Even if, initially, the couplings 
are taken to be constants, the renormalization process necessarily forces UiR to 
be momenta dependent all along FS. The counterterms, for this reason, form a 
continuum set in momenta space. 

Since our model is renormalizable, the counterterms originate naturally from 
the initial form of the hamiltonian. Here we neglect umklapp processes to start 
with. Then, the only divergent terms that arise in perturbation theory come from 
the interaction processes described by the U\r and U 2 r couplings as one can most 
easily verify. In view of this, in a first order approximation, we can disregard the 
perturbative terms coming from the other two couplings, namely U$r and U^r and 
their corresponding counterterms. As a result, the two sets of parallel patches of 
FS decouple from each other, and there is no ambiguity in locating a particle either 
at a solid line or at a dashed line patch instead. 

We begin by calculating first the electron self-energy in the vicinity of p± = kpR. 
Using conventional Feynman rules, the diagrams shown in fig. 3 produce 
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Figure 3: The diagrams used for the calculation of the renormalized self-energy. The 
ones with a cross represent the counterterm contributions, which regularize the diver- 
gences. Diagrams (d) and (e) cancel each other exactly. Diagram (f) produces the 
contributions directly associated with Z. 



where Um(p\\) = J dk\\Um{p\\, ku). The parameter A is the ultraviolet momen- 
tum cutoff with the corresponding energy cutoff given by Q = 2v FR X, and 2A is 
the length of each path along FS. In general, the renormalized coupling func- 
tions depend on three distinct momenta components parallel to the Fermi surface. 
However, in calculating the Hartree diagram, shown in fig. 3 (a), the vertex depend 
only on two different momenta and we naturally arrive at the definition of Um- 
In contrast, in the sunset diagrams of figs. 3(b) and 3(c), the renormalized vertices 
depend explicitly on three different momenta components along FS. In whatever 
way, in two-loop order, the tadpole diagram is in fact the only contribution from 
the self-energy which produce the renormalization of the Fermi surface. As a re- 
sult, for simplicity, in a zeroth order approximation, we take Um (p\\) = (2A)Um, 
neglect the momenta dependence of the vertices and rewrite the renormalized cou- 
pling in terms of the corresponding Urn's. For comparison, we show in Appendix 
1, how the diagrams and are modified if the full momenta dependence 
of the renormalized coupling functions are taken into account at all steps. 

In principle, since this a nonrelativistic system, there should be two scaling 
parameters in the problem: one for the energy (u) and another one for the mo- 
mentum (A) ^3]. By making an appropriate choice of the renormalization pre- 
scription, these two parameters do not mix with each other. Here we choose to 
work with the energy scale uj only. In doing this we implicitly assume that A — ► 
much faster than uj — > 0. Denning the renormalized one-particle irreducible func- 
tion T R (p±,p\\,p ;u), which is nothing but the inverse of the full single-particle 
Green's function, such that ReT R '(p± = kpR,p\\,Po = uj; uj ps 0) = uj, we find that 

(3A 2 - p 2 ) ~ ~ Q 
Z(p± = k FR ,pfu> « 0) = 1 - — " (U 2 1R + U 2 2R ) ln(-) + (4) 



and 



i*b = = z 1 (p ± = k FR ,pfutt o) (/^ + -^r ) + ( 5 ) 
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where we assumed, for simplicity, that the bare and renormalized masses differ 
only by a multiplicative factor, namely the charge renormalization function, i.e., 
m B = Z(p± = k FR ,p\\;uj « 0)m R . 



3 - Curvature and truncation in the renormalized 

FS 

Using perturbation theory, we calculate next the one-particle irreducible func- 
tions r(J (p, k, q;uj), which are essentially the renormalized two-particle interac- 
tion. Since our main intention in this work is to analyze the nature of the resulting 
renormalized FS, we restrict ourselves to scattering processes with two indepen- 
dent external momenta only. Employing appropriate renormalization group pre- 
scriptions such as 

r Su (Ph k bPo + ko = uj, k - p = uj;uj ~ 0) = -iU 1R (p\\,k\\;uj) , (6) 
T mn ip\\ > k \\ > Po + k o = w ~ 0) = -iU 2 r (p\\ , k\\ ; uj) , (7) 

and 

r 2i?TT(U) ( p \\ ' % Po + k ° = U] U ~ °) = °' ( 8 ) 

for p± = —k±_ = k FR , together with the perturbative expansions for the one- 
particle irreducible functions, we can relate the corresponding renormalized and 
bare coupling functions to each other. Following the same approximating scheme 
as before, with respect the dependence on the momentum component parallel to 
FS in the vertices of the various Feynman diagrams, and taking into account the 
RG conditions u)dU iB /duj = 0, we find the resulting RG equations 

1 d fj 2R (3A 2 - pfj „ _ (l7A 2 - 3p 2 ) _ 
P v™ = ^"af = Svr^ 2 A 2 UmUlR + 3847r^| D A 2 ^ + 



W -^I^1 P + ^9 P )+... 

FR^ ooi±/l u FR l 

(10) 



For comparison, we show in Appendix 2 how the RG equations for the corre- 
sponding renormalized coupling functions look like if we take into account the full 
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Figure 4: The flow diagram in the (U\ R: U2r) plane. 



momenta dependence at all stages. Finally, using eq.(5), and following a simi- 
lar procedure with the renormalized Fermi energy /i R (p\\,ui), we get the last RG 
equation of our interest 

_1 dfM R { 3A2 -P\\) f~ 2 ~ 2 \ ( SIU 1R \ 1 d (U m 
^ " Q W du ~ 1287^^0 V lR + Um ) \* R + 8^v FR ) &** W du [v FR 

(11) 

To determine the fixed points of the model, we calculate next the zeros of 
these RG equations for both renormalized couplings, and the renormalized Fermi 
energy. It then turns out that, aside from the usual infrared unstable Fermi liquid 
fixed point, and yet another nontrivial unstable fixed point, we get two infrared 
stable nontrivial fixed points which are, as we shall see, associated with non- Fermi 
liquid phases. For conciseness, we will only present the final expression for these 
IR stable fixed points 
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Figure 5: The anisotropic suppression of Z{p») in the present model as we approach 
the FS for three values of (ui/Cl). The black line is for (u;/f2) = 10 -5 , the light grey line 
is for (u/Q) = 10~ 7 , and the dotted line is for {ui/Q) = 10~ 9 . 

We observe that they depend upon py in a essential way. In view of that, the 
Fermi surface of the system also acquires a py-structure and deviates slightly from 
its initial flat form. This FS deformation comes out naturally from the renormal- 
ization process. This feature of the fixed points appears only in calculations up to 
two-loop order or beyond. In fig. 4, we show schematically the scaling trajectories 
in the (U 1R , U 2R ) plane. 

The nature of the electron liquid associated with the nontrivial fixed points 
can be inferred by the flow of the charge renormalization function Z(p\\,u). Since 
7 = (u/Z)(dZ/du), in the vicinity of the fixed point, we have that Z(p\\,u) scales 
as (uj/ST) 1 *^ with the anomalous dimension being simply 7 evaluated at those 
fixed point values. Note that it is always positive definite along FS. Calculating 
7*(pi|) explicitly, we find that, indeed, the charge renormalization function vanishes 
most rapidly at the center of the FS patch (fig. 5). This anisotropic suppression of 
Z(p\\,u) was also emphasized by Kishine and Yonemitsu within a Wilsonian RG 
approach [To] . 

Another particularly interesting feature is that, although the FS is IR stable 
within given boundaries, there is no guarantee that the physical FS is well-defined 
throughout the original patch. The nature of the resulting fermionic system can be 
determined by the behavior of the corresponding momentum distribution function 
n(p). If n(p) has an infinite slope at the renormalized FS, the state is metallic and 
resembles a Luttinger liquid. If it turns out otherwise that n(p) is perfectly smooth 
at kpR, the renormalized FS is truncated, and there appears a charge pseudogap 
at those points. This gapped state depicts an insulating spin liquid instead. To 
make the argument more quantitative, we use the momentum distribution function 
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Figure 6: The renormalized Fermi surface of the present model. The dashed line stands 
for truncation. 

calculated for a similar FS in ref. ^H]- In that work, it was shown explicitly that, 
for 1/2 < 7*(p||) < 1, the FS must be truncated at the corresponding p\\ values. 
Here this condition is fulfilled for |(pii/A)| < 0.41. After the elimination of the 
corresponding Fermi surface segments, the remains of the interacting FS are shown 
in fig. 6. This truncation scenario will be discussed in great length elsewhere |17j . 

4 - Conclusion 

In summary, we showed explicitly by a two- loop RG calculation that, even if the 
initial FS is entirely flat in two spatial dimensions, the true interacting Fermi 
surface becomes slightly curved as a result of interactions. Using field theoretical 
methods, we explored the existence of infrared stable nontrivial fixed points, which 
are associated with the non-Fermi liquid behavior in the low-energy sector. We 
showed that the fixed points vary continuously with the momentum p\\ along the 
Fermi surface. In addition, we argued that the criterion Z(u) — > as u — > does 
not suffice in determining the nature of the resulting non-Fermi liquid state. We 
called attention to a possible route of physical characterization of those states in 
terms of the momentum distribution function. For the pu values, in which the 
anomalous dimension "f(p\\) is such that 1/2 < 7*(pn) < 1, the renormalized FS 
suffers a truncation. This takes place, in our case, at the central region of the 
original FS patches. 

Not long ago, a new ARPES data for the high-temperature superconductors 
Bi2Sr2- x La x CuC>Q-s was reported showing some indication of spin-charge separa- 
tion and "Luttinger" liquid like behavior in the normal state of that compound [T5] . 
That data was, most recently, fitted consistently by a "Luttinger" liquid-like phe- 
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nomenology Those workers found that, contrary to what happens in one di- 
mension, the anomalous exponents vary strongly with momentum along the Fermi 
surface. This is in agreement with the findings of our work. We also note that, in 
the cuprates, there exists a pseudogap phase with a truncated FS around special 
k-values in momentum space. In view of that, it is therefore natural to ascribe to 
those gapped systems an insulating spin liquid nature instead [2T] . Again, it is 
suggestive to associate that result to our work. It is therefore important to explore 
this scenario further. In a even more general experimental framework, one could 
also try to apply pressure in a interacting metallic state to move the Fermi surface 
towards its critical condition inducing, in this way, a new kind of quantum phase 
transition [201 - This would, certainly, open more possibilities to test the limits of 
our results. 
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Cientifico e Tecnologico (CNPq) - Brazil. 

Appendix 1 

Using conventional Feynman rules, the contributions of the diagrams from figs. 3(b) 
and 3(c) for the self-energy and the non-interacting single-particle Green's func- 
tions associated with our model hamiltonian are 

^ {P\\,P± = k FR ) =-g^4 j D dk W d % (^J C/ lfl( fe ll>- fe ll+P||+?ll>?ll) 



*U 1R +P\\ +q\\) (po-v FR (p ± - k FR )) 



x 



and 



h| | A - vfr(p± ~ k FR ) -p -i5 \ k fn - v FR (p ± - k FR ) + p - id 



v FR (p± - k FR ) -po — iS J V v FR (p± - k FR ) +p -i5 



15) 



Sg c) (p hP± = k FR ) = -— L dk\\dq\\ U 2R {-k\\+p\\+q h k h q{) 



X U 2R (p\\,q\\,k\\) (po - v FR (p± - k FR )) 
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X 



n i A - vfr(p± - k FR ) - po - i5 ^ ^ /'Q - v FR (p ± - k FR ) +p -i5 



v F r(p± ~ k FR ) -po-iS 



v FR (p± - k FR ) +p -i5 



(16) 



where the domain of integration D3 is —A < fey < A, —A < q\\ < A, and 
—A < —k\\ + p\\ + q\\ < A. If, for simplicity, the coupling functions and the Fermi 
velocity are considered to be independent of the momenta components parallel to 
FS, our results follow immediately. 

Appendix 2 

Using perturbation theory together with our renormalization group prescriptions, 
and taking full account of the coupling functions dependence on the momenta com- 
ponents parallel to FS the RG equations for U 1R and U 2R , for general scattering 
processes, become 



CO 



9U lR (pi||,p 2 ||,P3||) 



du 



4vr 2 



dku 

D 1 \VFR 



UlR (p 1 \\,p 2 \\,k ll ) U 1R (pi|| +p 2 || - fc||,fc||,P3||) 



+U 2R (pi||,p 2 ||,A;||) U 2R {pi\\ +P2\\ - fc||,A;||,p3||) 



dk 

d 2 " \vfr 



UlR (Pl||,P3|| ~Pl\\ + k\\,P3\\) 



XUlR {k\\,P2\\,P3\\ ~Pl\\ + fc||) 



+ 64^4^ (Pl|hP2||,P3||) 



X 



E4 



'Pl||+P 2 ||,P3||+P4|| 



D 3 



dk\\dq\\ 



FR 



Uir +Pi\\ +q\\) 



x U lR {k\\,-k\\ +p 4 +q\\,q\\) + U 2R [p ih q\\, k\\) U 2R (-k\\ + Pi\\ + q\\,k\\,q\\) 



(17) 



and 



CO 



9U 2R {pi\\,P2\\,P3\\) 1 



du 



dku 



471-2 I JDi \ V FR 



Uir (pi||,p 2 ||,A;||) U 2R (pin +p 2 || - fc||,fc||,p 3 ||) 



+Ui R +p 2 || - fc||,%p3||) U 2R (pi||,p 2 ||,A;||) 



+ ^^2i? (P1||,P 2 ||,P3||) 
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4 f 

x ^2^pi\\+P2\\,ps\\+P4\\ / dk\\dq\\ . 2 

f=1 ./As \ V FR 



Uir (pi||,g||,-fc|| +Pi\\ +q\\) 



x U m (k h -k\\ +pi\\ +g||,g||) +U 2R {Pi\\,q\\,k\\) U 2R (-k\\ + pq +q\\,k h q ll ) 
where the domains Di, D 2 and D 3 are given by 

D 1 = 



(18) 



-A < k\\ < A 
-A < pin +p 2 || - h < A 



I>9 = 



-A < k\\ < A 

-A < p 3 || + fey < A 



-A < fc|| < A 
D 3 = ^ -A < q\\ < A 

-A < Pi || + q\\ - k\\ < A 

Once again if, for simplicity, we neglect the the dependence of the coupling func- 
tions with respect the components of the momenta along FS, consider scattering 
processes associated with only two independent external momenta and integrate 
over one of them, our simplified RG equations for the UiRS are readily reproduced. 
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